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Abstract 



In this paper we develop a Bethe approximation, based on the cluster vari- 
ation method, which is apt to study lattice models of branched polymers. 
We show that the method is extremely accurate in cases where exact results 
are known as, for instance, in the enumeration of spanning trees. Moreover, 
the expressions we obtain for the asymptotic number of spanning trees and 
lattice trees on a graph coincide with analogous expressions derived through 
different approaches. We study the phase diagram of lattice trees with nearest- 
neighbour attraction and branching energies. We find a collapse transition at 
a tricritical point, which separates an expanded phase from a compact phase. 
We compare our results for the transition in two and three dimensions with 
available numerical estimates. 

PACS numbers: 5.70.Fh, 36.20.Ey, 64.60. Cn 



The statistical properties of polymers have been one of the great challenges of statistical 
physics in the last decades (!],§]. A great understanding has been gained for the case of 
linear polymers, thanks to the simple but instructive description in terms of self-avoiding 
walks (SAWs) on a lattice, which is in turn amenable to many different treatments, ranging 
from Monte Carlo simulations to exact enumerations, from mapping to the 0(n = 0) spin 
model to field theoretical formulations ||. Less is known about the behaviour of branched 
polymers (BPs). In lattice statistical mechanics BPs can be efficiently modelled by lattice 
animals 0, i.e. by connected clusters of bonds. Analogously to linear polymers, BPs in 
solution can display a dense and a diluted phase, depending on temperature and quality of 
the solvent. The collapse transition occurs at the so called 9 point and in recent years 
has attracted a lot of attention, both for linear || and branched polymers |§. 

Recently, a variational approach based on the cluster variation method (CVM) has been 
introduced for linear polymers J7|, giving results in good agreement with the best numerical 
simulations in many situations: dense polymers (i.e. Hamiltonian Walks), diluted and 
self-interacting ones. The CVM |||9| is a closed form approximation, which finds wide 
applications in accurate investigations of the phase diagram of lattice spin systems [T(J . It 
is based on the minimisation of a variational free energy which is obtained by truncating 



the lattice cumulant expansion of the entropy ||. The largest clusters considered in the 
expansion determine the "order" of the approximation and are named maximal clusters. 
For example, standard mean-field and Bethe approximations are recovered by considering as 
maximal clusters respectively single sites and nearest-neighbours pairs. Most of the present 
knowledge about the phase behaviour of BPs has been gained through numerical approaches 
|| such as Monte Carlo simulations, exact enumerations or transfer matrix techniques (for 
an exception see, e.g., ref. (JTTJ ) . The possibility to extend to BPs the CVM scheme becomes 
therefore extremely appealing. 

In the present paper we introduce a CVM based variational approximation for lattice 
trees (LTs), that is lattice animals containing no closed loops. We show that the method per- 
forms very well when compared to situations where exact results are known as, for instance, 
in the case of the enumeration of spanning trees on a lattice. We then use it to explore the 
phase behaviour of LTs with branching energies and nearest-neighbour attractive interaction 
between non bonded monomers, in different energy and temperature regimes. 

Let us make first a rather important technical remark about our method. Being self- 
avoiding, LTs face strong non-local constraints on the lattice (neighbouring sites can be 
occupied by monomers that are extremely far apart on the tree). This is a serious problem 
when dealing with local approximations such as the CVM. In order to partially overcome 
this difficulty, we have introduced an intrinsic direction on each edge of a lattice tree, which 
can be visualised with an arrow. Each occupied site of the lattice is therefore characterised 
by a certain number of incoming edges (input edges) and by a certain number of outgoing 
edges (output edges). To ensure that the tree does not form loops we have required that 
(i) sites with no outputs occur with probability in the thermodynamic limit and (ii) all 
other occupied sites have exactly one and only one output. It is important to underline that 
conditions (i) and (ii) uniquely determine the orientation of the arrows on a given LT, so 
that the descriptions in terms of oriented LTs and non oriented LTs are equivalent. The 
advantage for our purposes is that we are able in this way to exclude configurations with 
an inner loop (which is a kind of long range constraint) by means of a local variable (the 
orientation of the arrow on each occupied bond). 

Within the CVM scheme, the Bethe approximation is obtained by choosing nearest 
neighbour (NN) pairs as maximal clusters. For a homogeneous and isotropic model with 
only single site and NN pair interaction terms, the variational (reduced) free energy density 
has the form || 

j( B ) _ _ |^Tr(p NN W NN ) + |Tr(p NN lnp NN ) - (q - l)Tr(p. In p.), (1) 

where q is the coordination number of the lattice, (3 = 1/fceT {ks = 1 from now on), Hnn 
is the contribution of a generic NN pair to the Hamiltonian and p, and pnn are the site 
and NN pair density matrices, respectively. For classical models the density matrices are 
diagonal and their diagonal elements are the probability of the corresponding configurations. 
The free energy must be minimised with respect to the density matrices, which must satisfy 
the condition of normalisation to 1 and compatibility (that is, p. must be obtainable from 
Pnn by a partial trace). 

The first problem we deal with are spanning trees (STs). A ST visits all sites of a 
lattice and therefore STs can be seen as a special subset) of LTs. In particular, 
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STs are believed to model compact branched polymers. For the moment we are interested 
in estimating the total number Afsr of STs on a lattice. To leading order, this number 
should scale with the number of sites N as Af ST ~ fi^ T , where fi ST is the so called connective 
constant. As the entropy per site s ST is given by s ST = In /i ST and there is no energy term, it 
follows that the reduced free energy per site f ST is simply related to the connective constant 
by fsr = — In jjLsT- The value of [x ST for a given graph is a non trivial number, which can 
anyway be calculated exactly. It is therefore an excellent starting point to test the accuracy 
of our method. 

The configurations of a site and a NN pair are classified according to the number of edges 
attached to each site and they are reported schematically in fig. [I]. In the following we will 
denote with the probability variable of a site configuration with i occupied edges (e.g. for 
a square lattice i takes values between 1 and 4) and with the probability variable assigned 
to a disconnected pair having % edges on one site and j edges on the other site. Similarly we 
will denote with /y the probability variable of a connected pair with an intrinsic direction 
from the site with % edges to the site with j edges. The free energy (|l|) takes then the form 

q f 9-1 9-1 9-1 - 1 f q ) 

fsr = E m e (i,j)e ij liae ij + 2j2J2 m f( i ^)fij ln fij ?-(?-!) \ zl m s{ i ) s i ln s< > (2) 

A {i,j=l i=0 j=l J U=l J 

where m s (i), m e (i,j) and mf(i,j) stand respectively for the multiplicity of site, disconnected 
pair and connected pair configurations (see fig. |I] for the actual values of m s (i), m e (i,j) and 
rrif(i,j)). The factor 2 in front of the contribution from the connected pairs configurations 
arises from the degeneracy associated with the direction of the arrow. In writing equation 
(0), moreover, we have implicitly assumed translational invariance in the system, a con- 
dition which in the thermodynamic limit is fulfilled. The normalisation and compatibility 
conditions on the density matrices can be written respectively as 

q-l q-lq-l ^ 

m e(hj)eij + 2j2J2 m f( i >j)fa = 1 ( 3 ) 

i,j=0 i=0 j=l 

and 

i- 1 . ,, - . , 

i = l,...q- 1 

(4) 




Conditions (Q) assures that single site and pair probability variables are defined consistently, 
so that by summing over all allowed configurations of one site of a pair, one should recover 
the probability s, associated to a single site configuration. Our problem is then to find 
the minimum of the CVM free energy f(p subject to the above constraints. In general this 
problem can be easily dealt with numerically with the help of an algorithm named numerical 
iteration method ||12||. In this special case of spanning trees enumeration however one can 
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first guess on the basis of numerical results and then verify by direct substitution that the 

(q — 2) q ~ i q 

solution takes the analytical form Sj = —, — — -, = ^Si-s?, fa = Sj + is, + i, which 

q(q — l)"? -1 q — 2 

gives for the reduced free energy 

= (g-2)ln[ g ( g -2)]-2( g -l)ln(g-l) 



from which one can derive the connective constant /x^? = exp(-f^))- 

It is interesting to remark that the same result for the reduced free energy (Q) can be 
obtained by exploiting the well-known relation between STs and the the Q- state Potts model 
]T3| . Indeed by applying the same approximation scheme to the Potts model and then taking 
the Q — > limit of the reduced free energy density calculated in f3 = ln(l + Q a ), where (3 
denotes the Potts coupling divided by ksT and < a < 1, one recovers expression (|5|). 

Actually, the exact number of STs on a <i-dimensional hypercubic lattice can be computed 
exactly JL4J]. The entropy is the logarithm of this number and, in the thermodynamic limit 



iV — > oo, the entropy per site is 



(2TT) d 

The large d expansion of (EJ) is 



\ p2it p2ir ( d \ 

SsT = TnZXd. J Q dkl '- J dk d \n \2d - 2j2cos(ki) j (6) 



SsT = M2(i )-l-^ + o(i) (7) 
Expanding (|) for large q we get 

4 s ; = \nq- — --L + o(\) (8) 
ST H 2q 2q 2 \q 3 J K ' 

With q = 2d the first three terms (the logarithm, the vanishing constant and the first power 
of 1/d) coincide. 

In fig. |2| we show a comparison between and the exact result s ST , which indicates 
that the accuracy of the approximation is rather good. This check is of utmost importance. 
Indeed, the CVM is an approximation on the entropy estimate of the system. Having an 
extremely good approximation of the entropy gives confidence in more complicated situations 
where energies are introduced into the model. The probability Pi that a randomly chosen 
site is connected to i nearest-neighbour sites can also be calculated exactly for STs on the 



square lattice @. The result is p 1 w 0.29454, p 2 « 0.44699, p 3 w 0.22239 and p 4 « 0.03608, 
which agree quite well with the values we obtain within our approximation (p^ = m s (i)si), 
i.e. px w 0.29630, p 2 w 0.44444, p 3 w 0.22222 and p 4 w 0.03704. 

As a further test, we have considered the case of directed spanning trees (DSTs) in d 
dimensions, where the exact result is known and it is simple. In this case a preferred direction 
along a lattice diagonal is chosen and all edges must have a positive component along this 
direction. This means that the output at each site is restricted to d possible directions. 
Moreover the output direction in the bulk of the system can be chosen independently one site 
from the other, as by construction loops cannot be formed. The number of DSTs is therefore 
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Most ~ d N (neglecting boundary terms). We have studied DSTs within our framework, 
selecting among all the single site and pair configurations for STs only those allowed by the 
directedness constraint (for instance, the number of single site configurations with % edges 

is m 8 (i) — d ^ ^ 1 ) ^ ) ^ r ^ s ' as com P are d to m s (i) — for STs). Similarly to the 

ST case, we have then written a variational free energy from ([[]). Minimising it numerically 
(subject to proper constraints), we have indeed verified that the exact result is recovered, 
i.e. fipsj, = d. 

In order to test the consistency of the scheme, we have studied the case of spanning trees 
with energies depending on the number of branchings at each site. To this end we have 
assigned a reduced energy penalty Ei to each site with i branches. In the limit Ex — > oo 
(and Ei = for i > 1) tips are not allowed, so that one should recover the case of space 
filling self-avoiding walk, i.e. Hamiltonian walks (HWs). We have verified that in fact when 
Ei — > oo the only configuration which have a non zero weight are those of a linear polymer. 
The resulting entropy moreover coincides with that of HWs in the Bethe approximation [|7| . 
Our Bethe approximation suggests a smooth crossover from STs to HWs, with no sign of 
discontinuity in the free energy. We have also verified, in the case q — 4, that the HW 
limit can be reached by sending simultaneously E 3 and E4 to infinity, thus forbidding the 
occurrence of branchings. 

We now turn to the problem of LT collapse. Just as for linear polymers, the collapse 
transition is driven by an attractive interaction f3 between nearest neighbour contacts^. The 
partition function of the so called t-model reads 

z N = j2t N (cy c (9) 

c>0 

where t N (c) denotes the number of trees with N sites and c contacts. Introducing a monomer 
fugacity z, the grand canonical partition function reads therefore 

00 

2=EE^W (io) 

JV=1 c>0 

where the first sum is over all possible number of sites in the tree. We then proceed anal- 
ogously to the ST case and write a pair approximation for the free energy of the system 
as 

I q 1- 1 

fix = - In z 2^771,(1)^ + -@ m e (i,j)eij + 

i=l 1 i,j=l 

q f 9-1 q ~ lq ~ l - 1 f q \ 

- { m e (i,j)e# Iney +^22J2 m f(^ j)fij ln fij ( ~ (<? ~ x ) m s (z)s; In s; \ (11) 
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i,j=0 i=0j=l I U=0 



where the symbols Sj, ey, have similar meanings as in eq. (^) (note that contrary to 
eq. (0) the possibility of an empty site must be considered in eq. ([□]))• The stable phase at 



A contact is defined as a pair of nearest neighbour vertices of the tree which are not linked by 
an edge 
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given (3 and z is obtained by minimising the free energy f£p subject to normalisation and 
consistency conditions on the probability variables (see previous discussion for STs). In this 
case it is not possible to express the solution in a simple, analytical form and one has to 
fully resort to numerical methods, e.g. the natural iteration method 

We report the complete phase diagram for d = 3 in fig. [| as a function of z and (3. There 
are two distinct phases: a zero density phase, where the average number of edges in a tree 
is finite, and a finite density phase, where instead this number is infinite. These two phases 
are separated by a transition line z c {(3) which could either be first order (corresponding to 
a finite jump in density) or second order (corresponding to a continuous change in density). 
The tricritical point where the two lines merge is the 9 point, /3g, i.e. the point where the 
collapse transition of the LT occurs. The structure of the infinite LT along the line z c ((3) is 
expanded for j3 < j3g and compact for {3 > /3g. We obtain (3g « 0.406 in d = 2 and j3g « 0.224 
in d — 3. Recent estimates from extensive Monte Carlo simulations on the collapse of lattice 
trees @ yield the values (3 e = 0.699 ± 0.052 and (3 e = 0.346 ± 0.017, respectively in d = 2 



and d = 3. 

The phase diagram for LTs appears to be identical to the one for SAWs (see e.g. fig. 2 

-j . In particular the 

numerical values j3e of the collapse transition of the polymers coincide and, in both cases, 
z c does not depend on /3, as long as f3 < /3g. This similarity is an intriguing result for 
which at the moment we don't have any plausible explanation. We have also investigated a 
generalisation of the partition function @ by including an energy term Ex which penalises 
configurations with tips. At any given E-y the phase diagram in (3 and z shares similar 
features with the phase diagram for LTs or SAWs. The difference is just restricted to a scale 
factor in z which is determined by the value of E\. This means in particular that the value 
of j3g does not depend on E\, at least at this level of the approximation. Similarly to the ST 
case, we have also verified that in the limit E\ —>■ oo, where tips are not allowed, the phase 
diagram of a linear polymer is recovered. 

It is interesting to remark that in the special case (3 = we numerically obtain for 
the connective constant of a LT \i LT = z c (0)~ l a result which agrees with the analytical 
expression 

(a- nta- 1 ) 

u ( s ) = Si J. do) 

^ LT (q - 2)(9" 2 ) 1 ' 



This expression was indeed derived in |T7[ by studying LTs on a Bethe lattice with coordi- 
nation number q and specializes to =6.75 for q = 4 and fi^J = 12.21 for q = 6. The 
most accurate numerical estimates [18] for the connective constant of LTs 5.14 
for the square lattice and jj, LT ~ 10.50 for the cubic lattice. 

In summary, we have introduced a novel variational technique to investigate the be- 
haviour of branched polymers. It is based on the cluster variation method with nearest 
neighbour pairs as maximal clusters. It can therefore be viewed as the natural formulation 
of the Bethe approximation for these systems and it provides the starting point on which 
systematic improvements can be made by applying the CVM recipe ||||| to clusters larger 
than the NN pair (the first step being the plaquette). Also, it should be possible to deal with 
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different branched polymer representations, such as weakly and strongly embedded lattice 
animals. 

The method yields exact results in the trivial case of directed spanning trees and ex- 
tremely accurate ones for the enumeration of spanning trees, where the exact result is still 
available. Introducing branching energies in the spanning tree problem we have shown that 
the Hamiltonian walk limit is correctly recovered. Finally, we have studied the collapse tran- 
sition of branched polymers, obtaining estimates for the point which compare reasonably 
well with the most accurate simulations. 

We thank A. Maritan for useful discussions. SL acknowledges EPSRC for financial 
support through a postdoctoral fellowship. 
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FIG. 1. Schematic representation of independent (a) site and (b) pair configurations, in the 
case of spanning trees. The continuous line represents the spanning tree; q = 2d is the coordination 
number of the lattice. The configurations drawn in the picture correspond to (a) S3, (b) e2,3 and 
h,2- 
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FIG. 2. Entropy per site as a function of coordination number q in the case of spanning trees. 
The continuous line corresponds to the Bethe approximation, the circles are exact results. 
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FIG. 3. Phase diagram of self-interacting lattice trees as a function of (5 and z for d = 3. 
The average number of bonds in a tree is finite (infinite) in region I (II). The continuous (dashed) 
line is a second (first) order transition. The cross marks the tricritical point (j3 e ~ 0.224 and 
z g 0.08192). 
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